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Introduction. 

Nonlinear differential equations and its systems are used to describe vari- 
ous processes in physics, biology, economics and so on. In nineteenth century 
it was shown that nonlinear equations are significant as much as linear ones, 
therefore there was the great progress of relevant mathematical theories, in- 
cluding the theory of nonlinear differential equations. 

Analytical solutions allow to obtain important patterns of model relation- 
ships, features of solutions behavior at different initial and border conditions. 
Exact solutions allow to test numerical calculations by one of the most reli- 
able ways. 

There are a lot of methods to look for exact solutions of nonlinear differ- 
ential equations: the inverse scattering transform [HE], Hirota method [3JH], 
the Backlund transform, the truncated Painleve expansion [5HH]. Here we 
present the generalization of modified simplest equation method, that was 
proposed to look for exact solutions of ordinary differential equations P|[TU]. 
Note, that one of advantages of this method is possibility of its complete 
automatization. 

The structure of the paper is as follows. In section 1 the modification of 
simplest equation method to look for exact solutions of nonlinear partial dif- 
ferential equations is presented. Using this method we obtain exact solutions 
of generalized Korteweg-de Vries equation with cubic source (section 2) and 
exact solutions of third-order Kudryashov-Sinelshchikov equation describing 
nonlinear waves in liquids with gas bubbles (section 3). 
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1 The simplest equation method to look for 
exact solutions of PDEs. 



The simplest equation method to look for exact solutions of ordinary 
differential equations is based on two main ideas: 1) on the usage of simplest 
equation with known solutions and 2) on the analysis of possible singularity 
orders of solutions |^I9HT3]. In papers [T3HT6] the simplest equation method 
is used to look for traveling-wave solutions of partial differential equations. 

Here we present the method analogous to modified simplest equation 
method, which allows to obtain the exact solutions of nonlinear PDEs that 
do not belong to the traveling-wave solutions class. 

We are looking for exact solutions of nonlinear PDEs of polynomial form 

P[u,u t ,...,4 k \u x ,...u®] = (1) 

We suggest that function u(x,t) can be presented via the solution of 
jjsimplest^ equation Z. As jjsimplest^ equation we use linear differential 
equations 

Zf = hZi l -V + ... + b^ 2 Z x + b { _ x Z (2) 
Z t = Cl Zt 1] + ■■■ + (H- 2 Z X + d-xZ (3) 

Sequentially the modified simplest equation method to look for exact 
solutions of nonlinear partial equation can be described in the following way. 

Step 1. Determination of possible orders of poles and zeros of traveling- 
wave solutions of equation studied. The singularity orders can be determined 
by step-by-step examination of leading terms or by power geometry methods 

wm. 

Step 2. Solution of PDE represents as function of jjsimplest equation^ 
solution taking into account orders of poles and zeroes, that have been de- 
termined at step 1 

u(x,t) = Y ( a i f i (Z) (4) 

i 

For example, if traveling-wave solution allows first-order pole, we can use the 
form 

u(x,t) = a + ai-^ (5) 
For n-order zero it's rational to use transformation 

u(x, t) = a + a%Z + . . . + a n Z n (6) 

and so on. 
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Step 3. Substituting into the initial equation and taking into ac- 
count relations (J2J) - (j3J) and setting coefficient at different powers of 
Z(x,t), . . . , Z(x, t) 1 ' 1 equal to zero, we obtain the algebraic system. Its so- 
lution determine coefficients czj and parameters of simplest equations bj, Cj 
(I = 1 . . ./ — 1). If it is necessary, the restrictions at parameters of initial 
equation can be done also. 

Step 4- Solution of equation fl2]) determines subject to restrictions, ob- 
tained at step 3, within I arbitrary functions of t, which by-turn are to be 
defined from equation ([3]) . Solution of original PDE expresses via solution of 
system ©-Q according to transformation (j3J). 

Note, this method can be completely automated by means of symbolic 
computations environment Maple. 

2 Exact solutions of generalized Korteweg-de 
Vries equation with cubic source. 

As a first example let us consider the generalized Kortewegde Vries equa- 
tion with cubic source 

u t + 3u 2 u x + 3(uu x ) x + u 3 + u xxx +p 1 u(u 2 +p 2 u+p 3 ) +p±u x +p 5 u xx +p 6 uu x = 

(7) 

The travelling-wave solutions of this equation allow first-order pole, so we 
are looking for exact solutions in the form 

u(x,t) — ax — , z = z(x,t) (8) 

z 

where 

Zxxx = hzxx + b 2 u x + b 3 u (9) 
Zt = C\Z XX + c 2 z x (10) 

Substituting transformation flS]) into original equation and using relations 
(fl9l- f|T0l) ). we obtain the algebraic system. Its solution gives us 1) restriction 
on coefficient of original equation 

p 6 =Pi + 2p 5 (11) 

and 2) parameters of relations flHl fTUj) : 

ai = 1, c 2 = pip 2 +P3 -P4, ci = Pi +P2 -P5, bi = -P2, b 2 = ~P3, h = (12) 

The type of solutions depends on sign p\ — 4p 3 , namely 
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at p\ - Ap 3 > 

C1A1 exp (Xix) + C 2 \ 2 exp (A2X) 



u(x, t) 



C exp ((pi + p 2 - Ps)p3t) + C\ exp (Aix) + C 2 exp (-A 2 x) 

X = X + (p 2 p 5 -Pl+P3~ P4)t, Ai )2 



2 .. w x . -P2 ± V 7 ^- 4 ^ 



at p 3 =p^/4 

Ad + C 2 (l + Ax) 



ti(x, £) 



Co exp (Ax) + Ci + C2X 



x = x + (p 2 p 5 - Pa - 3A 2 )t, x = x — (A (pi + p 5 + A) + p 4 )t, A = - — 
at - 4p 3 < 

(A 2 Ci - AiC 2 ) sin (Aix) + (Aid + A 2 C 2 ) cos (Aix) 



u(x, t) 



C exp (-A 2 x + (pi + p 2 - P5)p 3 t) + Ci sin(Aix) + C 2 cos (Aix) 
x = x + {p 2 p 5 -p 2 +p 3 - P4)t, Ai = - A 2 = — 

Here Ci, C2 and C3 are the arbitrary constants, non-zero one of which can 
be canceled. 

3 Exact solutions of third-order equation de- 
scribing nonlinear waves in liquids with gas 
bubbles. 

As the next example we consider the Kudryashov-Sinelshchikov equation 
for describing the pressure waves in liquid with gas bubbles 



u t + auu x + u xxx - (uu xx + uu x ) x - /3u x u xx - au xx = (13) 
Substituting transformations 

(14) 

u(x,t) = a + z n , z=z(x,t) (15) 
z xx = hz x + b 2 z, z t = c\z x + c 2 z (16) 
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into the original equation (|T5|) and equating coefficients at different powers 
of z and z x , we obtain that the modified simplest equation method results in 
exact solutions of equation ( fl3]) . if one of following conditions holds: 

P = -1, n — 2, bi — — 1, £> 2 = — , a — — - — , c x — — a - 1 - a, c 2 



4 2 ' * 4 

(17) 

2 , , a(2 + /3) 
a = -1, n = - - - , &! = -1, 6 2 = , a = 1, d = -a, c 2 = 

(18) 

« = -(TT^' n = 1,61 = ~TT^' = °' Cl = " (TT^i ^ = 

(19) 

a = —1, (3 = 2, n = 1/2, a = 1, b\ + &i — a = 0, 69 = 0, c\ = —a, c 2 = 

(20) 

Ot Ot ex 

p = 0, n = 1, 61 = -1, 6 2 = -, ci = --(1 + a ) - o - 1, c 2 = -(1 + tr) 

(21) 

Let us discuss resulting exact solutions at condition (TTTj) in detail. Func- 
tions u(x,t) are z(x,t) are to satisfy conditions 

u(x,t) = (a + 3)/2 + z(x,t) 2 
z t = (-a-l-a)z x + a{a + l)z/A (22) 
%xx z x -\- (XZ / 4 



So we obtain solutions of original equation (|T3|) 
• at P — — 2 a > -1 

m(x, t) = — h exp (zi) (Ci exp z 2 + C 2 exp {-z 2 )) 2 (23) 



at /3 = -2 a = -1 

I Q 

u{x,t) = ^— + e^(z 1 )(C 1 + C 2 (x-at)) 2 (24) 

at /3 = -2 a < -1 

m(x, t) = — h exp (21) [C\ sin z 2 + C 2 cos z 2 ) 2 (25) 
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where we use notation 
Z\ = ^1 + + a + - adj t — x, z 2 = — >J\(X + l\(at + t + crt — x) (26) 

Here C\ and Ci are the arbitrary constants. 

Solutions at conditions f[T8"j) . ( 1T^|) . f l2"Uj) or f[2"Tj) can be presented similarly. 
Note, that C2 = corresponds to the travelling-wave solutions. 



Conclusion. 

The modification of simplest equation method to look for exact solutions 
of PDEs is presented. The efficiency of this method was demonstrated by 
two different examples. Obtained solutions of the third-order Kudryashov- 
Sinelshchikov equation for describing the pressure waves in liquid with gas 
bubbles seem to be new. 
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